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Renormalization group analysis of the random first order transition 
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We consider the approach describing glass formation in liquids as a progressive trapping in an 
exponentially large number of metastable states. To go beyond the mean-field setting, we provide 
a real-space renormalization group (RG) analysis of the associated replica free-energy functional. 
The present approximation yields in finite dimensions an ideal glass transition similar to that found 
in mean field. However, we find that along the RG flow the properties associated with metastable 
glassy states, such as the configurational entropy, are only defined up to a characteristic length scale 
that diverges as one approaches the ideal glass transition. The critical exponents characterizing the 
vicinity of the transition are the usual ones associated with a first-order discontinuity fixed point. 
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In the ongoing search for a general theory of glass for- 
mation, the random first-order transition (RFOT) ap- 
proach has proven to be a strong candidate, establish- 
ing what appears to be an intricate mean-field (MF) de- 
scription of supercooled liquids and glasses. [lrl3j This MF 
treatment predicts a scenario with two critical temper- 
atures Td and Tk, the upper one Td being a dynami- 
cal singularity akin to the mode-coupling transition and 
the lower one Tk a thermodynamic ideal glass transi- 
tion characterized by a vanishing of the configurational 
entropy associated with the logarithm of the number 
of metastable states. This scenario has received sup- 
port from MF-like calculations on glassforming liquid 
models. HHJ The RFOT theory assumes that this MF 
picture retains some validity in finite-dimensional sys- 
tems and proceeds by accounting for ergodicity restoring 
and disappearance of metastability in supercooled liq- 
uids between Td and Tk through an entropy-driven nu- 
cleation process coupled with a mosaic view of the liquid 
configurations. [l|-|3|,|9j Testing the validity of this appeal- 
ing but still fragile scenario is of major interest. In addi- 
tion to computer simulations of model liquids, lfj,|ll| an- 
alytical work has so far been done in two directions, tak- 
ing the MF result as a starting point: instanton calcula- 
tions for the escape from glassy metastable states [12[|l3| 
and studies of disordered models with long-range inter- 
actions in the Kac limit. 1 12. Il4 | The latter in particu- 



lar have shown that both static and dynamic correlation 
lengths can be computed and that the MF infinite-range 
limit is not singular. Going beyond these approaches 
however requires a renormalization group (RG) treat- 
ment. 

We provide in this letter the first steps towards an RG 
treatment of glass formation beyond the RFOT MF the- 
ory. To this end, we consider the Migdal-Kadanoff (MK) 
real-space RG of a Ginzburg-Landau model which is com- 



monly taken to be in the "universality class" of struc- 
tural glass-formers as it generically displays the two- 
temperature scenario at the MF level. We use the replica 
formalism and, for convenience, we use hierarchical lat- 
tices, on which the MKRG is known to be exact. [151] 
Our starting point is the replica MF theory of structural 
glasses, in which one studies the distribution of putative 
metastable glassy states by introducing m — 1 copies (or 
replicas) of the same liquid system coupled with a small 
attractive interaction whose amplitude is set to zero af- 
ter taking the thermodynamic limit. 6, 16] By keeping 
the leading terms in the local order parameter, which is 
the similarity or "overlap" between different states, one 
obtains the following Ginzburg-Landau functional [131] 
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where q denotes the set of elements {q a b} (by con- 
struction, q aa — 0) and the overlap q a b(x) is physically 
associated with a local Debye- Waller factor character- 
izing molecular motion in the glass- forming liquid[13|; 
for simplicity, the only temperature dependence is 
taken in t — ~ ° , with T setting the temperature 
scale, while u,w,y > are considered as indepen- 
dent of temperature. This "real replica" method al- 
lows one to obtain the properties of the metastable 
states from the knowledge of the replica partition func- 
tion, Z(m) = fH ab T>q a b(x)exp(-S[q\). The mean 
free energy of a typical equilibrium state and the corre- 
sponding configurational entropy read respectively j3F — 



—d\ogZ(m)/dm and S c — — m 2 d (log Z(m)/m)/ dm. 
The number m of replicas should be analytically con- 
tinued to 1 in the equilibrium liquid phase and to a value 
less than one in the ideal glass phase, if present. [6| At the 
MF level, i.e. by looking for the uniform saddle-points 
of Eq. (flj , one finds that the order parameter q a b is zero 

2 

above a temperature Td, such that td = ^-, and that be- 
low Td appears another uniform solution with a replica 
symmetric (RS) structure q a b — q > for a ^ b. This 
solution, when plugged into the expression for F and S c , 
yields the properties of the metastable glassy states. Fi- 

9 2 

nally, at a temperature Tr- such that tx = 4f-, there is 
a RFOT with a coexistence between a zero-overlap phase 
and a high-overlap one, transition with zero latent heat 
and vanishing configurational entropy density. Below Tr , 
the system is in an ideal glass phase characterized by a 
nonzero overlap matrix and a value of m less than 1 . 
To go beyond MF, we consider a real-space MKRG ap- 
proach, which becomes exact on hierarchical diamond- 
like lattices, and apply it to a lattice version of the ef- 
fective Hamiltonian in Eq. ([T]). Such lattices are built 
iteratively by replacing each bond between sites by a 
fixed number of new bonds which, to mimic Euclidean 
d-dimensional lattices, is taken equal to 2 d . After n iter- 
ations, the volume of the system, which is equal to the 
total number of original bonds, is equal to 2 nd whereas 
the "distance" between the boundary sites is equal to 2 n 
bonds: this naturally fixes the length scale after n iter- 
ations as £ n — 2™. The procedure is illustrated in the 
inset of Fig. 1. The main advantage of this RG proce- 
dure is that the renormalized effective pair interaction 
between two sites at the n— th step of renormalization, 
W / n(q 1 J Q 2 )i satisfies a closed equation written in terms 
of the pair interaction W n _i(q 1 ,q 2 ): 

2- (d - 1) VF„(q\q 2 ) = 

log fl[dq ab exp {Wn-i(q\q) + V(q) + W„_i(q,q 2 )j 
a, b 

(3) 

where the labels 1 and 2 denote the value of two renor- 
malized sites from which emanate 2™( d_1 ) original bonds. 
At the nth iteration, the original lattice is replaced by a 
renormalized one where the unit length is £ n and the pair 
interaction between sites is W n (q 1 , q 2 ). A remaining ob- 
stacle is that even if the integration in Eq. ([3]) is purely lo- 
cal, it involves the components of a general mxm replica 
matrix with m continued to real values, and there is no 
general solution to this problem. However, the physics 
we aim at describing is related to that of first-order-like 
transitions; as a result, W n for generic q 1 and q 2 is ex- 
pected to grow rapidly with the number of iterations, as 
2 nd or 2™( <i ~ 1 ', so that when d > 2 the term in the expo- 
nential of the left-hand side of Eq. ([3]) becomes very large 
after a few iterations only. In consequence, we approxi- 
mate the full integral by a steepest-descent calculation. 
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FIG. 1: RG flow on the hierarchical lattice with d = 
3 for several initial conditions parametrized by T(= 
1.01,1.003,1.001,1.0001 from top to bottom): parametric 
plot of W n {q*,q*)/{2W n {q*,0)) versus 1/W n {q*,0). The two 
fixed points on the y-axis correspond to the RFOT (y = 0) 
and to the "normal" liquid (y — 1). Inset: Elementary step 
illustrating the RG on a hierachical lattice corresponding to 
d = 3. 



In addition, guided by the MF solution, we looked for 
p-step RSB saddle-points. Since we found that they all 
reduce to the RS saddle point, in the following we shall 
only focus on the latter. All these restrictions and limi- 
tations will be further discussed below. 
In the RS case, any matrix q is characterized by a 
single parameter q and, after introducing W n {qi,q 2 ) — 
to-1 an d V(q) — rrtri, the iteration equation sim- 
plifies to 
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and Tk, minimum which becomes the deepest one below 
Tk (the bare parameters are arbitrarily chosen such that 
tic = q(tK) = !)■ At the start of the RG flow, the "bare 
interaction" is taken equal to W (qi,q2) = (<Zi — q 2 ) 2 /2. 
It is easily derived from Eq. (j4]) that W n (0,0) does not 
flow and remains equal to zero (and so does the free en- 
ergyF„(0,0)). 

We first consider the liquid between Td and Tk- Asymp- 
totically, i.e. for large n, the system flows to a trivial 
disordered fixed point corresponding to a "normal" liq- 
uid and uncoupled replica. In this case, at each RG 
step, the minimum in Eq. Q is in q = 0, which leads for 
large enough n to W n (gi ,q 2 ) ~W n {q x , 0) + W n (q 2 , 0) for 
Qij'fe 7^ 0: boundary conditions then have no influence 
on the bulk. This behavior takes place along the RG 
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FIG. 2: Temperature dependence of the point-to-set correla- 
tion length £ps and of the penetration length £b in the liquid 
phase above Tk for a hierarchical lattice with d — 3. 
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FIG. 3: Renormalized configurational entropy density s c and 
surface tension v for a hierarchical lattice with d = 3. Main 
plot: Evolution with the RG scale £ n for T = T K + 0.5 12 . 
Inset: T-dependence of the values of s c and v just before the 
point-to-set scale. 



flow beyond a length scale £ps which depends on the ini- 
tial condition, here parametrized by T. More precisely, 
we define £ps = 2 nps from the number npg of itera- 
tions needed to make the value of q corresponding to 
the minimum in Eq. Q for both boundaries fixed in the 
metastable state q* drop below the value q* /2 (beyond 
this scale, the minimum is found in q = for all bound- 
ary conditions). From its definition, £ps corresponds to 
a "point-to-set" correlation length: }9Ul7| above £ps, the 
boundary conditions in the metastable glassy state do 
not affect the deep interior of the liquid whereas below 
£ps, they determine the state of the liquid. Exactly at 
Tk, one finds that the minimum sticks at q*^when the 
boundary conditions are fixed at q*: then, W n (q*,q*) 
stays equal to zero, just as W n (0,0). This corresponds 
to a first-order transition (a RFOT) with a coexistence 
between a liquid phase with q — and an ideal glass 
phase with q = <7*(Tr-). The fact that the value of Tk is 
itself not renormalized is a consequence of the minimiza- 
tion procedure (see also below). 

In Fig. 1, we illustrate the RG flow for several initial 
conditions in the case of a hierarchical lattice mimick- 
ing a 3-dimensional system (shown in Fig. 1). The be- 
havior is strongly reminiscent of that observed in a con- 
ventional first-order transition. [18| with the point-to-set 
length playing the role of the scale above which the low- 
T metastable phase disappears. We plot in Fig. 2 the 
T-dependence of £ps- It follows a power law, £ps{T) ~ 
(T — Tk) -1 - Interestingly, there appears to be another 
characteristic length scale £b, which we call the "pene- 
tration length" as it describes how far the "amorphous 
order" fixed by the metastable boundary condition pene- 
trates in the liquid; as such, £b seems to be related to the 
pattern repetition length introduced in [19J ] . Specifically, 
we compute it from the number of iterations rib required 



to make the value of q corresponding to the minimum in 
the RG equation for W n (q*,0) drop below g*/4. Above 
this length, the minimum is found at q — 0. The depen- 
dence of £b on initial conditions is shown in Fig. 2: it goes 
as (T — Tk)^' 3 ■ While this dependence might be specific 
to the structure of the hierarchical lattice, it is notewor- 
thy that the exponent 1/d is the same one obtained in 
[19j for the pattern repetition length and also coincides 
with the standard exponent for the "persistence length" 
near a first-order discontinuity fixed point. [20( Note that 
the asymptotic behavior of tpg and £b close to Tk can be 
also derived analytically (details will be presented else- 
where) . 

In order to obtain a description of the renormalized liq- 
uid on the scale £ n , it is useful to define and compute 
at the nth RG iteration both the configurational en- 
tropy, S c {£ ni T) = 0[F n (q*,q*)-F n {O,O)], and the in- 
terface free-energy T{£ n ,T) = F n (q*,0) - [F n (q*,q*) - 
F n (0,Q)]/2 between the liquid phase with q — and the 
glass phase with q = q* . Even though close to Tk S c 
is very small on microscopic scales, it increases by RG 
transformations as £^ whereas T grows as £^ _1 only. In 
consequence, S c and T eventually become of the same 
order. This happens on a scale of order £ps- At this 
point, the free-energy gain of having q — onsite becomes 
overwhelming compared to the free-energy cost due to a 
mismatch between the overlap values; ; the renormal- 
ized value of Upsilon then drops to zero, see Fig. 3, and 
the MF quantity Sc is no longer well-defined. [9| Note 
that on the lengthscale just below £ps the configurational 
entropy density s c = S c / '£„ linearly approaches zero as 
T — Tk, whereas the surface tension v = T /£^~ l remains 
nonzero at the transition, see Fig. 3.[25( The picture we 
obtain is that after nps RG iterations, the renormalized 
system is like a liquid at its "onset temperature" where 



the (renormalized) PS length is equal to one. The dif- 
ference with a normal liquid are the values of T and S c , 
which are very large compared to T, i.e. the liquid is 
at very low T compared to the typical scale of the inter- 
action. This-naively-suggests that the relaxation time 
could be obtained by assuming an Arrhenius law at the 
scale £ PS : logr ex T/T ex 1/(T - T K ) d - 1 . 
We now briefly discuss the situation below Tk- The cal- 
culation is similar to that performed above Tk except 
that when describing the equilibrium ideal glass phase 
with q = q*(T), the replica parameter m should also be 
optimized, which leads to m eq (T) < 1,[6] and that at 
the MF level, it is the liquid phase with q = which is 
metastable. We again find that the p-step RSB saddle- 
points reduce to the RS one (since m < L the extrem- 
ization on q now becomes a maximization 6] ) . We define 
a point-to-set correlation length by studying when the 
minimum involved in the iteration equation for W n (0, 0) 
increases above q* /2. We find that the equilibrium point- 
to-set correlation length defined for m = m eq (T) diverges 
as one approaches Tk from below with the same expo- 
nent as from above, i.e. as (Tk — T) _1 . 
From the above results, it appears that the physics in the 
vicinity of a RFOT to an ideal glass is controlled by a 
first-order discontinuity fixed point with standard expo- 
nents. Although the resulting physical picture is similar 
to the phenomenological one put forward in [l|,|9(, the val- 
ues of the exponents differs from [l|. Recent computer 
simulations of an atomic glass- forming liquid model. [Ill] 
find that some exponents are indeed standard, but others 
are not. It is therefore worth discussing the limitations of 
the present RG analysis. The starting replica Ginzburg- 
Landau model and the choice of hierachical lattices could 
of course be criticized, but the main potential shortcom- 
ing of the present study is the steepest-descent approxi- 
mation. We have argued that the latter is justified after a 
few iteration steps because the factor in the Boltzmann 
weight grows rapidly. In particular, if one starts with 
an arbitrarily strong random first-order transition at the 
MF level (this may not correspond to any actual glass- 
forming liquid but could nonetheless provide an interest- 
ing limiting case) and consider an initial condition deep 
in the glass phase, it is unlikely that the latter and the as- 
sociated random first-order transition would be destroyed 
by the few first iterations before minimization becomes 
justified. However, the parameters of the theory, such as 
the value of Tk , would clearly be renormalized by these 
first steps, so that minimization would start with effec- 
tive parameters in place of the bare ones.[26j Further- 
more, fluctuations may actually play an important role 
in determining the exponents of £b and consequently £ps 
since these involve the flow of S c that grows as (T— Tk)£„ 
and remains small until the scale £b- Another question- 
able aspect of our minimization scheme is the choice of 
saddle-points among RS and p-step RSB solutions only. 
Could we be missing relevant minima ? A way to get 



around the above problems would be to somewhat invert 
the replica treatment and infer from a replicated free- 
energy functional with a strong RFOT at the MF level 
an effective model with quenched disorder that could be 
directly studied by the MKRG on hierarchical lattices. 
This is a promising but arduous route, as there are at 
present no known disordered models displaying remains 
of the two-temperature scenario in finite dimensions. 22] 
(See however the recent proposal of a hierarchical ran- 
dom energy model in [23] and numerical simulations on 
a Kac version of the random orthogonal model 24] .) This 
nonetheless appears as the only means to check if, in a 
full-blown RG analysis, the trivial RFOT discontinuity 
fixed point found here is confirmed, replaced by a non- 
trivial one, or else if the RFOT itself is completely sup- 
pressed. 
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